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ABSTRACT 

Quantum  gates  and  simple  quantum  algorithms  can  be  designed  utilizing  the  diffraction  phenomena  of  a  photon 
within  a  multiplexed  holographic  element.  The  quantum  eigenstates  we  use  are  the  photon’s  linear  momentum 
(LM)  as  measured  by  the  number  of  waves  of  tilt  across  the  aperture.  Two  properties  of  quantum  computing 
within  the  circuit  model  make  this  approach  attractive.  First,  any  conditional  measurement  can  be  commuted  in 
time  with  any  unitary  quantum  gate  -  the  timeless  nature  of  quantum  computing.  Second,  photon  entanglement 
can  be  encoded  as  a  superposition  state  of  a  single  photon  in  a  higher-dimensional  state  space  afforded  by  LM.  Our 
theoretical  and  numerical  results  indicate  that  OptiGrate’s  photo-thermal  refractive  (PTR)  glass  is  an  enabling 
technology.  We  will  review  our  previous  design  of  a  quantum  projection  operator  and  give  credence  to  this 
approach  on  a  representative  quantum  gate  grounded  on  coupled-mode  theory  and  numerical  simulations,  all  with 
parameters  consistent  with  PTR  glass.  We  discuss  the  strengths  (high  efficiencies,  robustness  to  environment) 
and  limitations  (scalability,  crosstalk)  of  this  technology.  While  not  scalable,  the  utility  and  robustness  of  such 
optical  elements  for  broader  quantum  information  processing  applications  can  be  substantial. 

Keywords:  Linear  Optical  Quantum  Computing,  Volume  Holography,  Quantum  Gates 

1.  PHOTONIC  QUANTUM  ALGORITHMS  WITHIN  A  VOLUME  HOLOGRAM 

It  is  the  objective  of  this  manuscript  to  argue  that  the  construction  of  a  small,  lightweight,  field  deployable,  inex¬ 
pensive,  low-dimensional,  quantum  computer  is  feasible  using  existing  photonics  technology.  It  is  well  recognized 
that  photons  offer  great  promise  in  quantum  information  processing  (QIP)  given  their  robustness  to  decoherence. 
However,  it  is  this  very  resiliency  that  hinders  their  utility  in  quantum  computing.  Their  weak  coupling  with 
atomic  structures  leads  to  substantial  inefficiencies.  While  they  may  not  find  themselves  as  the  central  facet  of 
a  quantum  CPU,  they  can,  nevertheless,  be  an  integrable  part  of  the  quantum  information  processing  done  by  a 
quantum  computer  or  in  a  quantum  communication  system  where  repetitive  low-dimensional,  fixed  and  robust 
quantum  algorithmic  tasks  are  required  (e.g.  quantum  error  correction,  quantum  relays,  quantum  memory/CPU 
buss,  ...). 

Various  approaches  have  been  proposed  for  quantum  computing.  Perhaps  the  most  familiar  of  these  is  the 
quantum  circuit  model  (QCM).1  Here  the  challenge  is  to  identify  a  suitable  subset  of  a  universal  set  of  quantum 
gates  that  can  faithfully  represent  a  family  of  unitary  operations  on  d  qubits  (d-partite),  such  that  the  the 
number  of  gates  required  scale  polynomially  in  d,  whilst  the  dimension  grows  exponentially  as  2d  .  After  the 
unitary  evolution  the  output  needs  to  be  projected  onto  the  computational  basis  and  irreversibly  measured.  One 
realization  of  the  QCM  has  been  linear  optical  quantum  computing  (LOQC)  and  one-way  quantum  computing 
(OWQC),  which  are  also  known  as  cluster-state  quantum  computing. 2~5  The  essential  feature  of  each  of  these 
approaches  involves  either  a  non-linear  measuring  process  or  a  preparation  of  hyper-entangled  input  states  (|/A)), 
or  both.  Here  a  sequence  of  measurements  are  made  to  project  the  output  state  (| OUT))  onto  one  or  another 
of  the  computational  basis  states,  i.e  onto  a  mutually  unbiased  basis  (MUB).  In  this  manuscript  we  will  explore 
the  utility  of  phototnic  quantum  circuit  models. 

In  this  manuscript  we  suggest  that  in  principle  any  quantum  algorithm  can  be  encoded  within  a  single 
hologram,  and  that  in  practice,  it  has  been  found  that  one  requires  a  1  mm/N  thickness  of  holographic  material, 
where  ( N )  is  the  dimension  of  the  state  space.  This  thickness  will  ensure  high  fidelity  unitary  transformations. 

For  further  information  send  correspondence  to  W.A.M:  E-mail:  wam@fau.edu. 
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Nevertheless,  there  are  commercially  available  holographic  media,  e.g.  OptiGrate’s  photo-thermal  refractive 
(PTR)  glass,  that  can  support  faithful  holograms  with  a  thickness  approaching  a  few  centimeters,  thereby 
extending  our  considerations  to  10  to  20  dimensional  quantum  state  spaces.7  Every  application  of  LOQC 
to  quantum  gates  that  we  are  aware  of  requires  a  cascade  of  interferometers  which  is  highly  unstable  and 
impractical.  The  approach  here  will  “lock”  these  interferometers  within  a  tempered  piece  of  glass  that  is  resistant 
to  environmental  factors. 

In  this  paper  we  will  demonstrate  that  quantum  gates  and  simple  quantum  algorithms  can  be  efficiently 
simulated  by  the  diffraction  phenomena  of  a  photon  within  a  multiplexed  holographic  element.  In  Sec.  2  we 
briefly  outline  the  structure  and  function  of  a  volume  holographic  grating  that  will  be  referred  to  throughout 
this  manuscript.  In  Sec.  3  we  describe  the  photonic  quantum  eigenstates  we  use.  In  particular  we  use  the  space 
spanned  by  the  photon’s  linear  momentum  (LM)  as  measured  by  the  number  of  waves  of  tilt  across  the  aperture, 
i.e.  the  transverse  component  of  the  wave  vector  parallel  to  the  face  of  the  hologram.  In  Sec.  4  we  will  briefly 
review  our  previous  design  of  a  quantum  projection  operator.  In  Sec.  5  we  discuss  the  two  properties  of  quantum 
computing  within  the  QCM  that  make  this  approach  attractive.  First,  any  conditional  measurement  can  be 
commuted  in  time  with  any  unitary  quantum  gate  -  the  “timeless  nature  of  quantum  computing.”  Second,  photon 
entanglement  can  be  encoded  as  a  superposition  state  of  a  single  photon  in  the  higher-  dimensional  state  spaces 
afforded  by  LM.  Our  theoretical  and  numerical  results  indicate  that  OptiGrate’s  photo-thermal  refractive  (PTR) 
glass  is  an  enabling  technology.  Additionally,  we  give  credence  to  the  projection  operator  approach  of  Sec.  4 
by  applying  it  to  a  quantum  teleportation  (QT)  algorithm.  We  choose  this  algorithm  because  it  exhibits  both 
quantum  entanglement  and  the  timeless  nature  of  quantum  computing.1,8  We  used  coupled-mode  theory  and 
supportive  numerical  simulations  in  our  research,  all  with  parameters  consistent  with  PTR  glass.8,9  In  Sec.  6 
we  describe  the  construction  of  the  familiar  CNOT  gate  without  the  need  for  multiplexing.  Our  alternative 
approach  will  ordinarily  require  stacking  more  holographic  gratings.  In  conclusion,  (Sec.  7)  we  will  discuss  the 
strengths  (high  efficiencies,  robustness  to  environment)  and  limitations  (scalability,  crosstalk)  of  this  technology. 
The  holographic  approach  presented  here  is  not  scalable,  and  with  existing  technology  it  is  not  reprogrammable; 
nevertheless,  the  utility  and  robustness  of  such  optical  elements  for  broader  quantum  information  processing 
applications  can  be  substantial.  There  should  be  many  facets  of  a  quantum  computer  other  than  the  quantum 
CPU  that  could  take  advantage  of  this  technology,  including  error  correction,  QKD  relays,  quantum  memory 
management  tasks  and  the  like. 

2.  VOLUME  HOLOGRAPHIC  GRATINGS 

Volume  holography  is  used  today  for  2D  image  storage  utilizing  394 pixels / /im2 ,  which  consumes  only  1%  of  the 
theoretical  volumetric  storage  density  (1/A3).10  This  held  was  first  introduced  by  Dennis  Gabor  in  1948  and 
it  has  been  a  popular  research  held  ever  since  the  development  of  the  laser  in  1960.  Thicker  holograms  have 
more  precise  angular  selectivity,  i.e.  its  ability  to  differentiate  the  difference  between  two  plane  waves  separated 
by  a  small  angle,  and  under  certain  well-known  conditions  can  achieve  near  perfect  efficiency.11  A  hologram  is 
considered  a  volume  hologram  if  its  thickness  d  A2/A,  where  A  is  the  characteristic  period  of  the  index  of 
refraction  of  the  grating,  and  A  is  the  wavelength  of  the  light.  It  is  important  for  our  purposes  to  emphasize 
that  volume  holography  enables  higher  storage  densities,  and  under  suitable  recording  configurations  can  achieve 
near  perfect  efficiencies.  This  is  not  possible  with  thin  holograms  and  planar  gratings  (spatial  light  modulators) 
where  efficiencies  can  be  at  most  33%. 

Briefly,  the  transmission  volume  holograms  we  consider  here  are  formed  when  a  “signal”  wave,  (r\S)  = 
A{r)el^^r'>  is  directed  into  a  holographic  material  and  made  to  coherently  interfere  with  an  oblique  ’’reference” 
plane  wave,  \R)  as  illustrated  in  the  left  diagram  of  Fig.  1  for  N  =  3.  In  the  figure  the  “signal”  wave  is  a 
superposition  of  N  plane  waves, 

N  N 

<r1S)  =  ^(r15i)  =  ^e“‘e*'%  (1) 

i— t  i= 1 

where  cq  are  pure  phase  angles.  In  this  paper,  we  only  consider  planar  reference  waves,  and  the  signal  state  as  the 
superposition  of  plane  waves.  Ordinarily  the  signal  waves  will  have  variable  phase  and  amplitude  modulations. 
After  the  hologram  is  developed,  and  if  we  direct  the  identical  signal  wave,  (flS1),  into  the  hologram,  then  for  a 
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Figure  1.  The  left  diagram  shows  a  recording  of  a  volume  transmission  grating  by  the  coherent  superposition  of  a  plane 
reference  wave  |i?i)  and  a  linear  superposition  of  three  signal  waves  | S)  =  ezai\Si)  +  eI“2|/3S,2)  +  e*“3|S,3)  .  On  the  right 
we  show  the  function  of  the  hologram.  If  the  identically  oriented  signal  wave  IS)  is  sent  into  the  hologram  then  the 
reference  wave,  |-Ri),  will  be  reconstructed  in  the  diffraction.  The  diffraction  pattern  will  ordinarily  consist  of  higher 
order  diffracted  modes  parallel  to  the  signal  state.  However,  for  a  suitably  tuned  volume  hologram  perfect  efficiency  can 
be  achieved,  as  shown  in  the  right  diagram.8  This  is  why  we  constrained  the  signal  wave  components  to  a  cone  of  half 
angle  6S  centered  on  the  normal  to  the  hologram  face. 


perfectly  tuned  hologram,  the  reference  plane  wave,  (r|I?),  should  emerge  from  the  hologram  as  illustrated  in  the 
right  diagram  of  Fig.  1.  If  it  is  not  tuned,  then  other  diffracted  orders,  e.g.  modes  parallel  to  the  signal  states, 
may  emerge. 

How  does  one  form  a  perfectly  tuned  hologram  within  the  coupled  mode  theory?  We  have  shown  recently8 
that  near  perfect  efficiencies  can  be  obtained  if  (1)  the  hologram  thickness  is  tuned  to  its  optimal  thickness, 
(2)  if  the  each  of  the  signal’s  Fourier  wave  vectors  have  the  same  projection  onto  the  normal  to  the  hologram 
surface,  i.e.  they  all  lie  on  a  cone  with  half  angle  6S  as  shown  in  Fig.  2,  and  (3)  each  of  the  reference  plane 
waves  lie  on  their  own  distinct  cone  concentric  with  the  first,  with  half  angle  0r  and  centered  on  the  normal  to 
the  hologram  face.  We  will  also  consider  multiplex  holograms  wherein  multiple  independent  exposures  are  made 
within  the  holographic  material  before  it  is  developed.  We  demonstrated  using  coupled-mode  theory  that  if  the 
signal  waves  {5,i}i.i,2,...jv  form  an  orthogonal  set  under  the  L2  norm  in  the  plane  perpendicular  to  the  waves 
propagation  direction  (2),  i.e. 

{Si\Sj)  =  f  S*(x,y)Sj(x,y)dxdy  =  Sij,  (2) 

then  perfect  efficiency  can  be  achieved  for  each  of  the  signals.8 

A  volume  multiplexed  hologram  that  has  achieved  perfect  efficiency  (within  coupled-mode  theory12)  under 
the  “3+1”  conditions  outlined  above  provides  a  linear  map  between  signal  and  reference  modes.  Physically  it 
represents  a  projection  (or  redirection)  operator  or  signal  state  sorter,8 

V  =  |i?i}(Si|  +  \R2)(S2\  +  •  ■  •  +  \Rn)(SnI  (3) 

uniquely  identifying  each  pair  of  signal  and  reference  waves. 

Although  the  index  of  refraction  within  the  hologram  can  be  rather  complicated,  these  devices  are  strictly 
linear  optical  components.  Therefore,  the  diffraction  patterns  for  a  beam  of  photons  will  correspond  exactly  to 
the  probability  distribution  for  a  single  photon  in  the  beam.  For  the  remainder  of  this  manuscript  we  will  assume 
we  are  dealing  with  low  number  Fock  states.  We  will  describe  in  the  next  two  sections  (1)  how  we  encode  a 
d-partite  quantum  state  onto  a  single  photon  that  can  be  used  in  volume  holography,  and  (2)  how  the  unitary 
operation  representing  an  entire  quantum  algorithm  can  be  encoded  in  a  multiplexed  hologram. 
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3.  ENCODING  A  D-PARTITE  STATE  ON  A  SINGLE  PHOTON 


The  dimension,  N,  of  a  quantum  state  space  spanned  by  the  direct  product  of  d  qubits  grows  exponentially, 
N  =  2d .  Thus  to  represent  a  single  d-partite  state  by  a  single  photon  requires  vastly  more  quantum  numbers 
than  available  to  its  two  polarization  degrees  of  freedom.  However,  the  photon  can  be  characterized  by  its 
extrinsic  properties  as  characterized  by  its  wavefront  amplitude  and  phase.  The  potential  of  extending  photon- 
based  quantum  information  processing  and  quantum  computing  to  higher  dimensions  was  made  possible  in  1983 
when  Miller  and  Wheeler13  described  a  fundamental  quantum  experiment  utilizing  a  photon’s  orbital  angular 
momentum  (OAM),  and  when  Allen  et  al14  described  Laguerre-Gaussian  light  beams  that  possessed  a  quantized 
orbital  angular  momentum  (OAM)  of  lh  per  photon.  This  opened  up  an  arbitrarily  high  dimensional  quantum 
space  to  a  single  photon.13,15  Following  these  discoveries  Mair  et  al. 24,25  unequivocally  demonstrated  the 
quantum  nature  of  photon  OAM  by  showing  that  pairs  of  OAM  photons  can  be  entangled  using  the  non¬ 
linear  optical  process  of  parametric  down  conversion.  Shortly  thereafter,  Molina-Terriza  et  al.2'  introduced  a 
scheme  to  prepare  photons  in  multidimensional  vector  states  of  OAM  commencing  higher-dimensional  QIP,  with 
applications  to  quantum  key  distribution.  Recently  a  practical  method  has  been  demonstrated  to  produce  such 
MUB  states  using  computer-generated  holography  with  a  single  spatial  light  modulator  (SLM).16 

While  photons  with  specific  values  of  OAM  have  been  emphasized  recently  in  the  literature,  one  can  equally 
well  utilize  any  other  extrinsic  set  of  orthogonal  basis  functions  for  higher-dimensional  QIP,  e.g.  energy,  linear 
momentum,  angular  momentum.  While  OAM  states  respect  azimuthal  symmetry,  transverse  or  linear  momentum 
(LM)  states  respect  rectilinear  symmetry.  Independent  of  the  representation  we  use,  the  mutually  unbiased 
(MUB)  states  will  ordinarily  be  modulated  in  both  amplitude  and  phase.16 

While  the  advantage  of  higher-dimensional  QIP  lies  in  its  ability  to  increase  bandwidth  while  simultaneously 
tolerating  a  higher  bit  error  rate  (BER),17  such  transverse  photon  wave  functions  are  more  fragile  to  decoherence 
under  propagation  than  the  photon’s  spin  wave  function,18  and  the  divergence  of  the  states  in  propagation  may 
require  larger  apertures.  However,  since  the  spatial  distances  within  a  quantum  computer  are  small  this  is  not  a 
concern  here.  This  is  not  the  case  for  quantum  key  distribution.  Our  preliminary  results  show  that  the  efficiencies 
one  can  achieve  with  holographic  elements  exposed  by  LM  plane  waves  and  their  superpositions  (volume  Bragg 
gratings)  is  substantially  higher  than  that  for  photon  wavefronts  exhibiting  OAM.8,9  Multiplexed  volume  Bragg 
gratings  have  been  thoroughly  addressed  in  the  literature.'’12,19,20  Since  the  efficiencies  are  potentially  higher, 
and  volume  Bragg  gratings  are  well  understood,  we  therefore  will  restrict  our  attention  in  this  paper  to  the 
orthonormal  transverse  photon  states  of  LM. 

We  concentrate  in  this  section  on  the  physics  of  transverse  or  LM  states.  We  quantize  the  LM  within  the  2- 
dimensional  plane  of  the  face  of  the  hologram,  Fig.  2.  We  know  that  Bragg  gratings  can  obtain  perfect  efficiencies 
as  long  as  the  component  perpendicular  to  the  face  of  the  hologram  ( k± )  for  each  of  the  photons’  momentum 
are  the  same  i.e.  the  plane  wave  states  must  lie  on  a  cone.  In  addition  the  thickness  of  the  hologram  must  be 
carefully  chosen11,12  . 

If  we  are  to  represent  a  d-partite  state  by  a  single  LM  photon,  we  will  require  at  least  N  =  2d  LM  states 
around  the  cone  of  half  angle  9S  centered  on  the  normal  to  the  hologram’s  face.  As  a  result,  we  can  freely 
choose  as  our  computational  quantum  state  basis  any  2d  non-colinear  plane  waves.  In  this  case,  our  integer 
quantum  numbers  will  be  the  azimuthal  angles,  fa,  i  £  {1,2,  ...2d}  locating  the  wave  vector  of  the  photon 
around  the  cone.  Here  we  assume  the  hologram  surface  is  in  the  x-y  plane.  For  photons  of  wavenumber  k  these 
waves  correspond  to  a  transverse  linear  momentum  py  =  sin  (9s)Tik  =  Uk\ i,  with  components  pf  =  sin  (fa)Tik\\, 
Pi  =  cos  (fa)hk\\  and  p~  =  cos  (9)Tik  =  hk±  =  constant ;  respectively.  Here,  k\\  =  k(X/D)  is  the  magnitude 
of  the  transverse  component  of  the  wave  vector  of  a  plane  wave  of  wavelength  A  with  one  ’’wave  of  tilt”  (A /D) 
across  the  aperture  of  breadth  D. 

In  the  frame  of  the  hologram,  and  in  units  where  the  speed  of  light  is  unity,  the  components  of  the  3- 
momentum,  p  =  {px ,pv ,pz},  for  each  of  our  2d  photons  can  be  expressed  in  terms  of  their  wavenumber,  (k), 
i.e., 

Pi  =  hki  =  n  {kx,  k V,  kz}  =  hk  {sin  (9)  cos  (fa),  sin  (9)  sin  (fa),  cos  (9s)j  .  (4) 
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Figure  2.  These  Four  LM  states  (arrows)  that  lie  on  a  cone  are  used  to  generate  a  bipartite,  or  4-dimensional  signal  state 
space.  Here  we  chose  the  eigenstates  (each  a  plane  wave)  in  this  8-dimensional  space  with  eigenvectors,  | Si)  =  \9s,4>i),  i  £ 
1,  2, 3,  4,  here  (j>i  =  0,  (j>2  =  7r/2,  =  n  and  (f> 4  =  37t/2.  This  can  be  generalized  to  higher  dimensions,  2d,  by  introducing 

d  eigenstates  around  the  cone,  as  long  as  the  angular  selectivity  of  the  hologram  is  less  than  A(p  =  27r/2d.  The  wave 
vectors  are  defined  in  Eq.  4. 

These  N  =  2d  plane  wave  states,  (r|iSi)  =  elki'r  define  our  computational  basis  for  quantum  information  process¬ 
ing. 

MUB  =  {\S1),\S2),...,\Sn)}.  (5) 

Each  of  these  states  represents  a  transverse  Fourier  mode  of  a  photon;  they  are  orthogonal  under  the  L 2  norm 
((,S'j|;5'j)  =  5ij)  and  span  our  2d-dimensional  state  space. 

4.  REDIRECTING  STATES  IN  A  COMPUTATIONAL  BASIS 
FROM  REFERENCE  TO  STATE  SPACE 

In  Sec.  2  we  outlined  our  recent  progress  in  efficiently  sorting  a  single  photon  with  arbitrary  complex  wavefronts 
(Eq.  3). 8,9  And  we  showed  that  under  certain  conditions  perfect  efficiency)  within  coupled-mode  theory)  can 
be  achieved.  We  showed  that  a  multiplexed  hologram  could  function  as  a  d-partite  version  of  the  more  familiar 
qubit  polarizing  beam  splitter,  i.e.  a  single  optical  element  that  can  efficiently  sort  each  of  the  d-partite  states  in 
the  MUB ,  or  computational  basis  (Eq.  5).  As  we  argued,  thick  holographic  gratings  fortunately  produce  high 
diffraction  efficiency  in  the  first  order.11, 12, 19  If  several  predominant  diffracted  orders  are  required,  as  is  the  case 
for  sorting,  several  independent  fringe  structures  can  coexist  within  the  holographic  material.  Such  multiplexed 
holograms  have  been  used  for  multiple-beam  splitters  and  recombiners,19  and  more  recently  for  wide-angle  beam 
steering.21  In  a  recent  paper  we  described  such  a  MUB-state  sorter  based  on  a  multiplexed  thick  holographic 
element  constructed  from  commercially  available  photo  thermal  refractive  (PTR)  glass.8  Due  to  the  unique 
properties  of  PTR  glass  the  grating’s  thickness  can  approach  several  mm  and  be  highly  Bragg  selective.22  There 
is  evidence  that  such  sorters  can  be  highly  efficient,  with  Bragg  efficiency  >  95%. 21  Our  recent  simulations 
and  empirical  data  on  thick  Bragg  gratings  indicate  that  they  may  provide  an  adequate  solution  to  this  critical 
sorting  and  beam  steering  problem  for  “single  photon”  QIP. 

While  the  signal  state  sorter,  Eq.  3,  can  be  constructed  in  a  multiplexed  volume  hologram  under  the  strict 
conditions  outlined  in  Sec.  2,  we  could  just  as  well  construct  a  multiplexed  volume  Bragg  transmission  hologram, 
1Z,  that  recombines  a  beam.  i.e.  maps  the  reference  waves  onto  their  corresponding  signal  states.  Mathematically 
that  is  the  Hermitian  conjugate  of  Eq.  3, 

n  =  V^  =  |S'1)(R1|  +  IS2XR2I  +  ■  ■  ■  +  \Sn){Rn\.  (6) 
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This  can  be  constructed  physically  using  the  same  hologram  that  we  used  for  P,  but  instead  of  directing  the 
signal  waves  to  the  hologram  to  reconstruct  the  corresponding  reference  wave,  we  could  direct  a  reference  wave 
which  will  reconstruct  the  desired  superposition  of  signal  wave  states,  i.e.  the  signal  beam,  IS").  For  example, 
two  identical  multiplexed  holograms  ( V )  stacked  back-to-back  would  represent  the  identity  operator,  X,  on  this 
state  space.  Any  signal  wave  used  to  record  the  hologram  emerges  unchanged  from  the  holographic  sandwich; 
however,  within  the  holographic  material  it  will  undergo  a  complicated  evolution  into  the  reference  states, 
and  back.  Though,  this  identity  gate  may  not  be  particularly  useful  in  its  own  right,  it  does  help  illustrate  our 
point.  We  will  find  in  the  next  two  sections  that  these  holographic  state  redirection  elements  (Eq.  3)  can  be 
useful  in  the  construction  of  quantum  computing  gates. 

5.  QUANTUM  TELEPORTATION  IN  GLASS 

To  best  illustrate  how  a  quantum  algorithm  can  be  encoded  within  a  single  hologram  we  will  focus  our  attention 
on  quantum  teleportation  ( QT ).  This  three  qubit  gate  lives  in  an  8-dimensional  state  space.  It  exhibits  all  the 
essential  features  we  need  to  address,  i.e.  entanglement  and  conditional  measurements.  The  approach  here  will 
generalize,  in  principle,  to  any  quantum  circuit  model.  We  will  adhere  as  closely  as  we  can  to  the  notation  found 
in  Nielsen  and  Chuang.1 

The  standard  experimental  realization  of  QT  utilizes  a  pair  of  entangled  qubits,  |/?oo),  and  an  arbitrary  qu&it 
state  vector  to  be  teleported,  ip,  where, 


/loo)  = 

^=(|00>  +  |11» 

(7) 

ip  = 

a|0)  +  P\l) 

(8) 

The  QT  circuit  consists  of  a  CNOT  gate,  a  Hadamard  gate  as  well  as  two  conditional  measurements,  one  of 
which  is  fed  into  a  Pauli  X  gate,  the  other  into  a  Pauli  Z  gate  as  shown  in  Fig.  3.  We  exploit  the  ” timelessness” 
of  quantum  computing  by  using  the  principle  of  deferred  measurements  and  commute  the  two  measurements  in 
Fig.  3  with  the  quantum  X  and  Z  gates.  We  arrive  at  an  equivalent  circuit  model  for  QT  where  all  measurements 
are  made  in  the  end  as  illustrated  in  Fig.  4.  We  focus  in  this  section  on  this  equivalent  tripartite  QT  circuit. 
Essentially  QT  amounts  to  (1)  preforming  a  unitary  transformation  on  the  input  tripartite  state  followed  by,  (2) 
two  measurements, 

\OUT)  =  UQT  \IN),  (9) 

where 

I  IN)  =  \ip)  ®  |/30o>.  (10) 

This  unitary  transformation,  Uqt,  can  be  decomposed  into  a  product  of  four  sequential  unitary  transformations 
as  illustrated  by  the  dashed  box  shown  in  Fig.  4.  Evolution  proceeds  from  the  left  to  the  right.  The  input  state 
first  encounters  (1)  a  CNOT  gate  on  the  first  two  qubits  with  unitary  transformation,  Ucnot ,  followed  by  (2) 
a  Hadamard  transformation  on  the  first  qubit  with  unitary  operator,  IAh ,  followed  by  (3)  a  Pauli  X  gate  on  the 
third  qubit  conditioned  on  the  second  qubit  with  unitary  transformation,  Hex,  and  finally  (4)  a  Pauli  Z  gate 
on  the  third  qubit  conditioned  on  the  first  qubit  with  unitary  transformation,  Ucz-  We  can  therefore  write  the 
overall  QT  unitary  transformation  as  a  product  of  these  four  transformations, 

Uqt  =  Ucz  Ucx  Uh  Ucnot  (11) 

This  8-dimensional  state  space  can  be  interpreted  as  a  direct  product  of  3  qubit  states  (Eq.  10)  with  entan¬ 
glement  in  the  last  two  qubits  (Eq.  7).  However,  we  can  also  represent  the  tripartite  state  as  a  pure  state  in 
an  8-dimensional  state  space  without  loss  of  generalization.  In  this  way  the  entire  input  state  (Eq.  10)  can  be 
represented  by  a  single  photon  in  an  8-dimensional  LM  eigenspace  as  described  in  Sec.  4.  We  just  need  to  define 
a  correspondence  between  a  state  in  the  computational  basis  and  one  of  the  LM  eigenstates  (Eq.  5).  We  make 
this  identification  below. 
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Figure  3.  A  quantum  teleportation  circuit.1  Here  the  8-dimensional  tripartite  input  state  | IN)  =  \ip)  ®  |/3oo)  is  the  direct 
product  of  the  arbitrary  qubit  state,  \ip)  =  l/-\/2  (|0)  +  |1))  we  wish  to  teleport,  and  a  maximally  entangled  pair  of  qubit 
states  |/3oo)  =  1/72(1 00}  4-  1 1 1) ) .  This  circuit  is  composed  of  an  initial  CNOT  operation,  followed  by  a  Hadamard  gate, 
and  two  Pauli  transformations  (A'  and  Z)  conditioned  on  two  measurements. 


It  is  convenient  for  our  analysis,  and  our  application  to  holography,  to  represent  the  equivalent  QT  algorithm 
in  matrix  notation.  Here  we  choose  the  standard  representation  of  the  qubits, 


and  |1)  = 


The  two  input  states  can  then  be  expressed  as  vectors  in  a  2  and  4  dimensional  state  space;  respectively, 

a 

P 


IA)o)  — 


1 


(  1  \ 


72 


0 
0 

V  1  7 

Therefore  the  input  state  can  be  expressed  as  a  vector  in  an  8-dimensional  state  space, 

(  a  \ 

0 

0 

a 

p 

0 
0 

\p  J 


I  IN)  =  7)  ®  \Poo)  =  -j= 


(12) 


(13) 


(14) 


(15) 


Following  the  definitions  of  Sec.  3,  we  can  freely  choose  the  following  computational  basis  vectors  in  this 
8-dimensional  state  space  and  their  correspondence  to  each  of  eight  “signal”  LM  plane  waves: 


(r|<S'1)  =  eifcl'r"' 


/1\ 
0 
0 
0 
0 
0 
0 

V  0  J 


,  ( f\S2)=eik ** 


t  0  \ 
1 
0 
0 
0 
0 
0 

V  0  J 


(r\S$)=eik8'f! 


(  0  \ 
0 
0 
0 
0 
0 
0 

V  1  J 


(16) 
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Figure  4.  A  equivalent  quantum  teleportation  circuit.1  We  use  the  “principle  of  deferred  measurement”  and  commute 
the  time  ordering  of  the  conditional  measurements  with  the  last  two  Pauli  gates.  This  emphasizes  the  “timeless”  feature 
of  quantum  information  processing.  The  QT  circuit  thus  becomes  a  single  unitary  transformation  followed  by  two 
measurements  and  the  transported  state  (top  wire  — >  bottom  wire). 


The  quantum  numbers  of  the  “signal”  state, 

I  Si)  =  |fcf>  (17) 

represents  the  projection  of  the  linear  momentum  of  the  photon  onto  the  face  of  the  hologram,  i.e.  the  transverse 
components  of  the  LM.  There  would  ordinarily  be  two  components  of  this  momentum,  | Si)  =  .  /cj ) ,  hence 

two  quantum  numbers.  However,  here  each  of  these  eight  “signal”  plane  waves  lie  on  a  cone  of  half  angle  9S. 
Therefore,  the  first  quantum  number  is  constant  for  all  states, 

k ^  =  k^~  =  constant ,  Vi.  (18) 

We  define  the  transverse  quantum  numbers  by  equally  spacing  the  states  around  the  cone’s  perimeter  with  (f>\  =  0 
and  A <fi  =  7t/4.  This  is  a  generalization  to  eight  plane  waves  of  the  four  wave  configuration  illustrated  in  Fig.  2. 
Therefore,  one  can  say  each  state  is  described  by  two  quantum  numbers,  9  and  (f> ;  however,  we  hold  9  fixed 
so  that  our  holograms  can  achieve  perfect  efficiency  (within  coupled  mode  theory)  in  modeling  a  given  unitary 
transformation.  This  will  be  described  below. 

The  input  state  (Eq.  15)  is  then  a  linear  superposition  of  four  of  the  eight  computational  basis  vectors;  and 
correspondingly,  can  be  represented  by  a  linear  superposition  of  four  plane  wave  state  vectors, 

I  IN)  =  -j=  (a|Sr)  +  a|S4>  +  P\S5)  +  0\Ss)) .  (19) 

This  state  has  a  wavefront  with  both  amplitude  and  phase  variations.  In  a  similar  vein  we  can  define  a  corre¬ 
sponding  set  of  eight  “reference”  plane  waves  that  will  be  used  to  construct  our  holograms.  Here  we  define  them 
as  we  did  for  the  “signal”  waves,  on  a  seperate  concentric  cone  of  half  angle,  9r  ^  9 s, 

(r\Ri)=eiK'1"r,  <f1  R2)  =  eiK2r,  •••,  (r\R8)  =  eiK°'r .  (20) 

Similarly,  each  of  these  reference  waves  live  on  a  second  cone  of  half  angle  9r  ^  9S  and  are  equally  spaced,  but 
offset  with  respect  to  the  signal  waves,  around  its  perimeter  with  </> i  =  7r  and  A (j)  =  7t/4.  It  is  important  to  note 
that  the  exact  placement  of  the  eigenstates  around  the  cones  and  relative  relationships  between  the  reference 
waves  and  the  signal  waves  are  unimportant  within  the  couple-mode  theory.  Any  configuration  can  yield  perfect 
efficiency.  However,  in  practice  it  is  preferable  to  make  the  angular  separation  between  the  signal  and  reference 
waves  as  large  as  possible,  since  this  can  lead  to  higher  angular  selectivity  and  higher  efficiencies.  This  is  not  a 
concern  in  this  paper. 

We  can  use  these  basis  vectors  to  construct  each  of  the  four  unitary  transformations  shown  in  the  dashed 
boxes  of  Fig.  4,  thereby  yielding  a  representation  of  the  overall  QT  unitary  transformation,  Uqt,  given  by  Eq.  11. 
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In  particular  we  find, 
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as  well  as, 
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(21) 


(22) 


Therefore,  in  our  computational  basis  the  unitary  transformation  representing  the  QT  can  be  obtained  by 
substituting  Eqs.  21-22  into  Eq.  11, 


°QT  =  7I 


(  1  0 

0  1 

0  0 

0  0 

-1  0 

0  -1 

0  0 

\  0  0 


0  0  0 

0  0  0 

0  10 

10  1 
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0  0  0 
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-10  1 


0  1  0  \ 
0  0  1 
10  0 
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(23) 


This  unitary  transformation  of  the  equivalent  QT  circuit  shown  in  Fig.  4  accomplishes  the  desired  transformation, 


I  IN) 


top 

wire 


bottom 
two  wires 


\OUT)  =  Hqt  \IN) 


top  wire 


middle  wire 


^  ’ 
bottom 
wire 


(24) 


thereby  teleporting  the  state  function  from  the  top  wire  to  the  bottom  wire. 

We  describe  how  to  encode  the  unitary  transformation,  Hqt  into  a  single  holographic  element.  By  construc¬ 
tion,  the  rows  of  Hqt  form  an  orthonormal  basis  that  spans  the  8-dimensional  state  space.  Each  of  these  rows 
can  be  expressed  as  a  superposition  of  the  computational  basis  state  vectors  of  Eq.  16.  For  this  example,  each 
row  requires  a  superposition  of  exactly  two  of  the  eight  computational  basis  state  vectors.  Therefore,  each  row 
can  be  represented  as  a  linear  superposition  of  the  corresponding  (Eq.  16)  plane  wave  states  -  two  plane  waves 
each  for  this  particular  quantum  circuit.  For  example,  if  we  consider  the  state  vector  representing  each  of  the 

rows  of  the  matrix  in  Eq.  23,  we  can  construct  each  of  the  the  eight  QT  signal  states,  |£j)  that  we  will  use  in 

the  successive  recordings  in  the  hologram, 

|Sr)  =  (\Sr)  +  |Si»  <*=►  U^.  (25) 

|E2)  =  -±=  (|S8>  +  |<?2»  <*=►  Uqt2j  (26) 
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|S3)  =  72(I^6>  +  |54))  ^  UQT3,j 
\^)  =  ±(\S5)  +  \S3))  <*=►  U^QTij 

\^5)  =  ^(\s7)-\s1))  <*=►  uy6J 
|E6)  =  -±=  (|58)  -  |S2»  <*=►  u'QTf,. 
|S7)  = -^  (|S6)  -  |54))  UqT7j 

|E8)  =  -j=  (|S6)  -  |S3»  <*=►  Uqts^- 


(27) 

(28) 

(29) 

(30) 

(31) 

(32) 


By,  inspection  we  can  reconfirm  that  these  eight  basis  vectors  form  an  orthonormal  basis  spanning  the  8- 
dimensional  state  space. 

We  construct  a  coordinate  system  where  the  z-axis  is  orthogonal  to  the  face  (the  x-y  plane)  of  the  hologram.  A 
faithful  holographic  representation  of  Uqt  can  be  achieved  by  multiplexing  eight  recordings  within  the  holographic 
material.  The  hologram  is  recorded  so  that  each  row  of  the  unitary  matrix,  Uqt ,  is  used  to  generate  its  own 
volume  holographic  grating.  A  qubit  gate  such  as  the  Hadamard  gate  (H)  is  represented  by  a  2  x  2  matrix  and 
would  ordinarily  require  two  recordings.  A  2-qubit  gate  such  as  the  controlled-not  gate  CNOT  is  represented 
by  a  4  x  4  matrix  and  would  ordinarily  require  four  independent  recordings  (this  will  be  discussed  in  the  next 
section) .  The  qutrit  QT  circuit  is  represented  by  a  8  x  8  dimensional  matrix  (Eq.  23)  and  requires  the  following 
eight  sets  of  independent  recordings,  between  each  of  the  “signal”  states  in  Eqs.  25-32  and  its  corresponding 
“reference”  plane  wave  of  Eq.  20  as  follows: 


1.  The  first  recording  is  made  by  the  coherent  superposition  of  “signal  state  (r|Si)  and 

2.  The  second  recording  is  made  by  the  coherent  superposition  of  “signal  state  (r|E2)  and  (r|i?2); 


3.  The  eighth  and  final  recording  is  made  by  the  coherent  superposition  of  “signal  state  (r|Eg)  and  (r|i?g) 

The  recording  geometry  for  the  case  of  the  fifth  recording  is  illustrated  in  Fig.  5. 

Since  (1)  each  of  the  eight  signal  states  lie  on  a  single  cone  with  axis  normal  to  the  face  of  the  hologram  and 
form  an  orthonormal  basis  spanning  the  state  space,  and  (2)  each  of  the  eight  reference  states  lie  on  another 
concentric  cone  and  span  an  isomorphic  8-dimensional  state  space,  the  oct-tuple  multiplexed  hologram  (within 
coupled  mode  theory)  faithfully  represents  the  desired  transformation, 

Uqt  =  ^  (|i?i)<Si|  +  \Ri)(S7\  +  \R2)(S2\  +  \R2)(S8\  +  \R3)(Si\  +  \R3)(S6\  +  \R4)(S3\  +  \Ri)(S6\- 
|ik)<Si|  +  \R5)(S7 1  -  \R*)(S2\  +  \Re)(S8\  -  \Rr)(St\  +  \R7)(Se\  -  \Rs){S3\  +  \Rs){Ss |) . 

However,  this  operator  maps  quantum  states  in  ‘signal’  space  (inner  cone  of  Fig.  5)  onto  corresponding  states  in 
the  reference  space  (outer  cone  of  Fig.  5).  In  order  to  map  this  output  back  into  the  corresponding  states  in  the 
signal  we  would  need  to  take  the  diffracted  output  of  the  hologram  described  by  Eq.  33  as  input  to  a  “redirection” 
as  described  by  Eq.  6  in  Sec.  3.  This  can  be  constructed  in  a  similar  fashion  by  eight  recordings,  wherein  the  \th 
recording  is  made  by  the  coherent  superposition  of  the  \th  signal  plane  wave  basis  state  (rib)  =  elki'r  with  the 
1th  reference  plane  wave  basis  state,  (r\Ri)  =  elKi'r,  yielding, 

8 

RQT  =  ^2\Si)(Ri\.  (33) 

i= 1 
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Figure  5.  Volume  holographic  design  of  the  fifth  of  eight  independent  recordings  used  to  generate  the  QT  unitary  trans¬ 
formation,  Uqt,  in  PTR  glass. 


Stacking  the  redirection  hologram  behind  the  hologram  for  the  desired  transformation  will  yield  the  desired  QT 
operation. 

QT  =  TZqt  Uqt  =  ^  (|5i)(Si|  +  \Si){S7\  +  |S2)(S2|  +  |S2><S8|  +  |53}(54|  +  S3)(S6\  +  |S4)<S3|  +  |54)(S5|- 
|55)(5i|  +  \S5)(S7\  -  \S6)(S2\  +  \S6)(S8\  -  |S7)<S4|  +  \S7)(S6\  -  \S8)(S3\  +  |S8)<S5|) 


6.  THE  CNOT  GATE:  STACKING  INSTEAD  OF  MULTIPLEXING 

One  alternative  to  multiplexing  is  to  make  single  recordings  in  each  of  many  holograms  and  to  stack  the  holograms. 
Here  we  provide  a  second  illustrative  example  with  a  specific  design  of  a  quantum  CNOT  gate  compatible  with 
PTR  glass.  This  gate  is  realized  by  stacking  four  holograms,  and  we  will  describe  each  below.  The  CNOT  gate 
is  a  two  qubit  gate.  Therefore  the  dimension  of  the  state  space  is  4-dimensional.  While  this  state  space  can  be 
constructed  as  a  product  space  of  qubits  by  utilizing  the  polarization  states  of  two  correlated  photons,  it  can  also 
be  represented  by  a  single  LM  photon  in  a  4-dimensional  state  space.  The  CNOT  gate  can  be  constructed  with 
a  single  photon.  Following  the  arguments  of  the  previous  section,  we  freely  choose  four  independent  plane  waves 
lying  on  the  cone  shown  in  Fig.  6.  We  associate  to  these  independent  transverse  LM  modes  the  four  orthogonal 
quantum  state  vectors,  |5i),  |iS2),  |<S3)  and  |54).  Any  state  vector,  | ip),  in  this  4-dimensional  state  space  can  be 
written  as  a  linear  superposition  of  these  states, 


with, 


|V>)  =  a|S1}  +  /?|S2}  +  7|S3}  +  <5|S4}, 


(34) 


M2  +  |/3|2  +  |7|2+H2  =  l- 


(35) 
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S4 


Figure  6.  Volume  holographic  design  of  the  4-dimensional  CNOT  gate  in  PTR  glass.  The  gate  can  be  constructed  by  a 
stack  of  4  LM  gratings,  or  by  a  stack  of  two  multiplexed  gratings. 


If  we  let  the  z-axis  be  orthogonal  to  the  face  ( x-y  plane)  of  the  hologram,  the  four  volume  holographic  gratings 
are  recorded  by  a  suitable  superposition  of  the  set  of  four  signal  plane  waves, 

(r|5i)  =  exp  (iki  •  rj  ,  (r\b)  =  exp  (zfo  •  r) ,  (r^)  =  exp  ^ ik 3  ■  rj  ,  and  ( f\d )  =  exp  (ik^  ■  r) ,  (38) 

and  four  reference  waves, 

(r|i?i)  =  exp  (i«i  •  r) ,  (r|i?2)  =  exp  (in 2  •  r) ,  (r|i?3)  =  exp  (in 3  •  r) ,  and  (r|i?4)  =  exp  (1^4  •  r) ,  (39) 

as  shown  in  Fig.  6. 

The  hologram  is  recorded  so  that  each  row  of  the  unitary  matrix  of  the  CNOT  gate  is  used  to  generate  its 
own  volume  holographic  grating.  For  a  2-qubit  gate  such  as  the  CNOT  gate  we  would  ordinarily  require  four 
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recordings;  however,  since  the  first  two  bits  are  just  an  identity  matrix  we  need  only  two  layers  to  transform  the 
signal  states  into  the  desired  reference  states.  In  addition  to  one  holographic  recording  per  dimension  of  the  state 
space,  we  also  require  the  conjugate  of  the  each  grating  (two  in  the  case  of  the  CNOT  gate)  in  order  to  transform 
the  diffracted  reference  waves  from  the  reference  waves  back  into  the  desired  signal  states.  In  particular,  the 
CNOT-gate  constructed  from  four  holographic  gratings  stacked  together  as  is  shown  in  Fig.  6. 

1.  The  first  grating  is  recorded  with  the  two  coherent  plane  waves  corresponding  to  states  \S3)  and  |i?4)  . 

2.  The  second  grating  is  recorded  with  the  two  coherent  plane  waves  corresponding  to  states  £4)  and  \Ft.3)  . 

3.  The  third  grating  is  recorded  with  the  two  coherent  plane  waves  corresponding  to  states  |i?4)  and  IS4)  . 

4.  The  fourth  grating  is  recorded  with  the  two  coherent  plane  waves  corresponding  to  states  and  IS3)  . 

The  four  gates  will  not  diffract  the  first  two  signal  states  |£l)  or  | £2).  However,  the  first  two  gratings  redirect 
the  two  signal  states  5,3 )  and  (£4)  into  |i?4)  and  \R3)  respectively,  in  accordance  with  the  Pauli  X-gate, 

*-(!  »)■  (40) 

The  first  hologram  is  equivalent  to  the  operator, 

Hi  =  |5i)(5i|  +  |S2)(S2|  +  \Ri)(S3\  +  IS4XS4I,  (41) 

and  the  second  hologram  recorded  with  signal  plane  wave,  (f]£3)  and  reference  plane  wave,  (r|i?4),  is  equivalent 
to  the  operator, 

U2  HSrXSrl  +  \S2)(S2\  +  \R4)(R4\  +  \R3)(S4\.  (42) 

After  the  signal  state  ,  | IN),  passes  through  the  first  hologram  the  basis  vectors  change  from  {£1,  £2,  £3,  £4} 
to  the  orthonormal  basis,  {£1,  £2, 1?4,  £4}.  This  explains  why  the  third  term  in  Eq.  42  is  entirely  within  the 
reference  space. 

While  these  two  recordings  could  have  been  made  in  a  single  multiplexed  hologram,  we  recover  the  same 
function  by  stacking  the  two  together,  thereby  generating  the  CNOT  operation, 

Ucnot  =  U2UX  =  |£i)<Si|  +  |£2)(£2|  +  |-R3)(£4|  +  l^4)(£3|-  (43) 

However,  the  output  of  these  two  stacked  holograms  are  the  reference  states,  |i?2),  I-R3)  and  I-R4).  In  order 
to  redirect  these  back  to  the  proper  signal  states,  we  require  the  redirection  operator  similar  to  Eq.  6.  This  can 
be  accomplished  by  recording  a  third  hologram  with  the  states  \R3)  and  |£3).  The  third  hologram  is  equivalent 
to  the  operator, 

U3  =  |£iX£i|  +  |52)(£2|  +  \S3)(R3\  +  \R4)(Ra\-  (44) 

Similarily,  the  fourth  hologram  is  recorded  with  the  states  \Ra)  and  (£4)  and  the  fourth  hologram  is  equivalent 
to  operator, 

UA  =  |£i)<£i|  +  |£2)(£2|  +  |£3)(£3|  +  |£4)(iZ4|.  (45) 

Therefore,  the  combination  of  the  four  stacked  volume  holograms  have  the  desired  action  -  the  CNOT  gate. 

CNOT  =  U2-U3-U2-U1  =  |£i)(£i|  +  |£2>  <£2 1  +  |£4)(£3|  +  |£3)(£4|.  (46) 

One  can  apply  these  principles  to  design  a  universal  set  of  quantum  gates,  as  well  as  simple  quantum  algorithms 
such  as  QT. 

The  advantage  of  stacking  the  holograms  is  that  one  can  make  the  hologram  thicker,  thereby  increasing  the 
efficiency;  however,  achieving  and  maintaining  the  proper  alignment  should  be  more  problematic.  By  multiplex¬ 
ing,  we  would  need  two  holograms,  each  with  two  independent  recordings  in  them.  The  first  would  be  equivalent 
to  the  last  two  holograms  in  Fig.  6,  while  the  second  would  be  equivalent  to  the  first  two  and  would  just  redirect 
the  reference  beams  into  their  corresponding  signal  states.  The  first  two  recordings  are  complementary  to  the 
second  two  -  thus  in  some  sense  we  are  recording  the  “square  root”  of  the  CNOT  gate. 
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7.  DISCUSSION:  ADVANTAGES  AND  LIMITATIONS 


Constructing  simple  quantum  algorithms  and  quantum  gates  in  volume  holograms  provides  substantially  greater 
optical  stability  than  the  equivalent  optical  bench  realization.  For  LOQC,  this  stability  is  the  overarching 
advantage.  Often  quantum  operators,  e.g.  the  simple  projection  operator  given  by  Eq.  3,  require  a  cascade  of 
interferometers  where  the  output  of  one  interferometer  is  used  as  the  input  of  the  next.5, 23,26  Therefore,  as 
the  dimension  of  each  state  space  increases,  it  becomes  exceedingly  hard  to  stabilize  and  is  simply  impractical 
beyond  two  qubits.  Other  approaches,  such  as  crossed  thin  gratings  lack  the  efficiency  needed  for  QIP.  The  device 
proposed  here  can  potentially  achieve  this  in  a  single  piece  of  glass  without  the  problem  of  misalignment.  The 
technology  presented  here  can  potentially  replace  ’’fixed”  optical  components  for  a  broad  spectrum  of  classical 
and  quantum  photonics  experiments. 

The  primary  limitation  of  volume  holographic  QIP  is  that  it  is  not  scalable.  Experience  shows  that  mul¬ 
tiplexing  requires  approximately  1  mm  per  recording  of  the  state  space  to  achieve  high  fidelity,  and  in  QIP 
applications  this  scales  exponentially  with  the  number  of  qubits.  But  then  again,  we  are  not  aware  of  a  realistic 
scalable  quantum  computer  to  date.  Secondly,  the  holograms  discussed  here  are  write-once  holograms  and  can¬ 
not  be  erased.  Therefore  the  algorithm  is  “fixed”  into  the  holographic  element.  While  there  are  re-recordable 
holographic  media,  none  that  we  know  of  has  the  specifications  to  outperform  PTR  glass  for  the  applications 
discussed  in  this  manuscript.  This  is  hardly  a  prescription  for  a  quantum  CPU;  however,  as  mentioned  in  the 
paper,  this  technology  might  be  integral  to  complete  QIP  systems  where  smaller  d-partite  operations  are  needed 
on  a  routine  basis,  e.g.  a  quantum  memory  bus,  quantum  error  correction  circuit,  QKD  relay  system,  etc... 

While  we  have  extensively  analyzed  volume  holographic  gates  and  algorithms  using  coupled-mode  theory, 
paraxial  wave  equation  simulations  and  finite-difference  time  domain  simulations,  we  have  not  analyzed  the 
engineering  particulars  of  this  device.  For  example:  (1)  how  many  independent  writes  of  orthogonal  states  into 
a  holographic  element  can  be  made  in  the  PTR  glass  before  cross-talk  or  saturation  between  the  modes  becomes 
a  limiting  factor?  (2)  Is  it  difficult  to  stack  the  holograms  due  to  the  enhanced  angular  selectivity  of  the  volume 
holograms?  And  (3)  what  is  the  maximum  number  of  recordings  in  a  multiplexed  PTR  hologram  that  can 
be  reasonably  be  achieved?  In  this  sense,  we  are  well  along  in  understanding  these  devices  from  a  theoretical 
prospective;  however,  we  are  at  the  very  beginning  experimentally. 
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